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A new theory is presented for the development of circular polarization as radio waves propagate 
through the turbulent, birefringent interstellar medium. The fourth order moments of the wavefield 
are calculated and it is shown that unpolarized incident radiation develops a nonzero variance in 
circular polarization. A magnetized turbulent medium causes the Stokes parameters to scintillate 
in a non-identical manner. A specific model for this effect is developed for the case of density 
fluctuations in a uniform magnetic field. 



I. INTRODUCTION 



Circular polarization observed in the radio emission from pulsars and quasars is not understood. Although the 
emission mechanisms for these two classes of source are quite different, both have the common feature that the 
emission is due to highly relativistic particles in a magnetic field, for which the polarization should be predominantly 
linear with a circular component of order the inverse of the Lorentz factor of the radiating particles. As discussed 
further below, this intrinsic component of circular polarization does not account for the observations. We propose 
that the circular polarization is imposed as a propagation effect, due the scintillations in the interstellar medium £] ] 
which account well for many of the observed, and otherwise unexplained, time variations in pulsars and quasars 
The underlying physics for this alternative explanation is presented in this paper. We propose to discuss the details 
of the astrophysical applications elsewhere. 

Scintillations are attributed to scattering of the radio waves off density inhomogeneities associated with turbulence 
in the ISM. The data indicate a power-law model for the turbulence, with the power-law index consistent with the 
Kolmogorov value (f3 ~ 11/3 in the notation used here) [Q. The dataset on pulsars is sufficiently large to allow 
mapping of the turbulence across the Galaxy implying much stronger scattering at low Galactic latitudes, where 
most pulsars are observed, than at high Galactic latitudes, where most quasars are observed. A planar wavefront 
becomes rippled as it traverses the region of turbulence. Two length scales play a central role in the theory: the 
Fresnel scale, rp — (XD/2-k) 1 ^ 2 , where A is the wavelength of the radio wave and D is the distance between the 
observer and the scattering region, and the diffractive scale, rdiff, over which the fluctuations in the phase decorrelate 
due to the turbulence. Physically, rdiff characterizes the sizes of the ripples, and rp characterizes the size of the 
coherent patch an observer can see on the unrippled wavefront when only the geometric phase difference is taken into 
account. "Weak scattering" corresponds to rp <C rdiff , when an observer sees a single coherent patch that is slightly 
tilted (image displacement) and slightly convex (focussed) or concave (defocussed). "Strong scattering" corresponds 
to rp 3> r^ff, when an observer sees many coherent patches of size rdiff within an envelope of size r rc f = rp /rdiff, 
and multipath propagation occurs Intensity variations in strong scattering are induced by both diffractive effects, 
due to interference between the coherent patches, and refractive effects, caused by focussing (defocussing) of the ray 
bundle due to phase curvature across the scattering disk. There is a transition between strong scattering at lower 
frequencies and weak scattering at higher frequencies. For most pulsars the transition frequency is higher than the 
frequencies for which data are available, and for quasars the observations span the expected transition frequency, 
~7GHz§. 

Inclusion of the interstellar magnetic field implies that the ISM is birefringent, and propagation of radiation depends 
upon its polarization. In a homogeneous birefringent medium, radiation separates into the two oppositely polarized 
natural wave modes, with the phase difference between them increasing linearly with propagation distance. The two 
wavefronts corresponding to the two modes become systematically displaced from each other with increasing distance. 
The natural modes in the ISM are circularly polarized to an excellent approximation, and the birefringence results 
in Faraday rotation of the plane of linear polarization of any incident radiation. The amount of Faraday rotation is 
parametrized in terms of the rotation measure, RM, which is defined such that the phase difference between the two 
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modes is RMA 2 ||. Inhomogeneity combined with birefringence implies that the wavefronts are both rippled and 
displaced from each other. One implication is that when the wavefronts are recombined there is a random component 
in the phase difference, and the resulting stochastic Faraday rotation is characterized by a variance in RM ||. 

The main point made in the present paper is that scattering in the magnetized ISM necessarily leads to a component 
of circular polarization (CP) in the observed radiation. This arises because any lateral displacement of the wavcfront 
implies that the ripples are not superimposed when the two wavefronts are combined. As a result, there are alternative 
patches of excess right hand (RH) and excess left hand (LH) circular polarization on the image in the observer's plane. 
In this paper we present a detailed theory for such scintillation-induced CP. We argue that the predicted features of 
the CP are sufficiently promising to warrant the development of detailed models for the observed CP in pulsars and 
quasars based on it. 

CP is observed in both pulsars and in some compact extragalactic sources JTT] |L1|] . Most of the data on pulsars 
are integrated over many pulses, and the integrated pulse profile typically shows relatively small CP. However, there is 
evidence that in at least a few pulsars for which individual pulses can be studied, the CP is relatively large in individual 
pulses and varies from pulse to pulse such that the integrated value is much smaller than the typical value. There is no 
satisfactory explanation for the C P fl5|l . A small 0.1% to a few %) but significant degree of CP is observed in some 
compact extragalactic sources pT|-|l4[. The suggested interpretations include the intrinsic polarization associated with 
synchrotron radiation [fl6| , and partial conversion of linear into circular polarization due to ellipticity of the natural 
wave modes of the cold background plasma (l^] or of the relativistic electron gas itself ]l8|-p0| . However none of these 
suggested interpretations has proved satisfactory in accounting for (a) the frequency dependence, (b) the temporal 
variations, and (c) the magnitude of the observed circular polarization The explanation proposed in this paper is 
as outlined above. Specifically, for a source with zero intrinsic CP seen through a turbulent magnetized plasma, there 
is a variable CP which corresponds to a zero average of the Stokes parameter V , (V) = 0, but a nonzero variance, 
(V 2 ) 7^ 0. Our initial objective is to use the theory of scintillations in a magnetized plasma p2|-|2l| to calculate (V 2 ). 

The magnitude of the expected value of CP needs to be of the same order of magnitude as the observed CP for 
the theory proposed here to be relevant. For pulsars, the CP can be several tens of percent, but it may be that 
some of this CP results from birefringence in the pulsar magnetosphere itself, which we do not consider in detail 
here. For the most extreme case for quasars, the observed CP can be several percent, which is relatively high because 
there is independent evidence suggesting that the varying (scintillating in our interpretation) part of the source is 
only a modest fraction of the entire source. Hence, the observations, in the most extreme cases, suggest that the CP 
of the scintillating component can be as high as a few tens of percent. In the theory developed here, most of the 
terms that contribute to the CP are very small, of order the ratio of the cyclotron frequency in the ISM to the radio 
frequency (typically ~ 10 -8 ), and can be of no practical interest. However, the effect on which we concentrate can give 
arbitrarily high CP. This effect is due to birefringent refraction causing an angular separation between the emerging 
rays in the LH and RH polarized wave modes. The displacement of the centroids of the LH and RH images increases 
linearly with the distance from the screen where the birefringent refraction occurs. It is possible for the LH and RH 
images not to overlap, resulting in patches of 100% CP. The angular deviation required to produce relatively large 
CP is determined by the ratio of the characteristic size of the scintillation pattern divided by the distance between 
the observer and the screen where the birefringent refraction occurs. Although the angular separation between the 
rays is always extremely small, we argue elsewhere |26| that observed gradients in RM imply birefringent refraction 
through a sufficiently large angle to satisfy the criterion that observable CP be produced. 

Our specific assumptions are explained in §2, where the wave equation is reduced to the form used in treating 
the scattering. The mutual coherence of a polarized wavefield is derived in §3, and is shown to reproduce some 
known results [^7],^8| . In §4 the second-order correlations of the Stokes parameters are derived from the fourth-order 
moment of the wavefield, and explicit solutions are obtained in the thin-screen approximation. In §5 we discuss 
scintillation-induced CP. The conclusions are presented in §6. 



II. PROPAGATION THROUGH A MAGNETIZED PLASMA 



In this section the propagation of radiation through a magnetized stochastic medium is related to its effect upon 
the electric field of the wave. We start with propagation through a homogeneous weakly anisotropic medium and then 
generalize to include the effect of inhomogeneities in the scattering medium. For a weakly inhomogeneous medium 
the two wave modes are assumed to be transverse to a zeroth approximation (the isotropic limit) and the degeneracy 
between the two transverse states of polarization is broken by the weak anisotropy, to a first approximation. In the 
first approximation the two modes have slightly different refractive indices, and this approximation suffices to treat 
Faraday rotation and all the effects of interest here. 



The wave equation projected onto the transverse plane is 
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k) = 0, 



(2.1) 



where A a (ui, k) is the wave amplitude, K a "(oj, k) is the dielectric tensor and the greek indices run over the two 
transverse co-ordinates. We write 



(2.2) 



where the angular brackets denote the mean and 8K a @ denotes a fluctuating part with a mean of zero, (SK a P) = 0. 
The transverse components of the dielectric tensor may be expressed in terms of the Pauli matrices, and for the 
average part we write 



(K af3 ) = 

where the sum over A = [/, Q, U, V] is implied, with 
= S = (o 1 ' a Q = 0-1 
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(2.3) 
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In the discussion here we neglect any dissipation, which implies that we retain only the Hermitian part of K a @ , so 
the Ka are all real. 

The two modes are labelled a = ±. The mode a has wavenumber k a = n a uj/c, where n a is the refractive index for 
waves in mode a, and polarization vector e CT . It is straightforward to solve for the dispersion relations, which are 



kl = \{K I + a[Kl + Kl + K*} 1/2 } 



(2.5) 



The quantity Kj is the dielectric constant in the isotropic approximation, with Kq + 
anisotropy limit. The polarization vectors in the two-dimensional transverse plane are 



(K Q + a \K% + Kfj +K^] 1/2 ,Ku + iK 



err = 



{2[Kl + Kl + Kl\ 1/2 (K Q + a[Kl + Kl + K* v ] ^) } 



2l l/2 u l/2- 



Ky <C Kj in the weak 



(2.6) 



One is always free to orient the coordinate axes such that Kjj — 0. On doing so, equation (2.6) may be written in the 
form 



(l,*T«r), 



T„ = 



Kx 



Kq + a \Kq + K v ] 



1/2' 



T+T- = -1, 



(2.7) 



where T a is the axial ratio of the polarization ellipse in the mode a. In the case of Faraday rotation, the wave modes 
are circularly polarized, which corresponds to \T a \ = 1. The circular polarizations are 



e rA = -J=(l,±l), 



(2.8) 



where r and I refer to right- and left-hand respectively. In this paper we consider only the circularly polarized 
approximation explicitly, but the general theory is valid for any T a . 

In scattering theory it is conventional to make the parabolic approximation to the wave equation p9j , |30] ] . This 
approximation is related to the paraxial approximation in geometric optics, in the sense that there is a favored ray 
direction (the z-axis in our case) and that only small deviations from it are considered. In the parabolic approximation 
the wavefield is written as the product of a fast- varying term e lh ° z and a term u a (z, r) that is assumed to be a slowly 
varying function of z in the sense that its second derivative with respect to z may be neglected. The difference 
between fc_|_ and k_ is introduced explicitly by writing k a = k + aSk, which defines 8k. In a homogeneous medium 
the two phase factors ik a z include the mean phase ikz, corresponding to the average over the two modes, and the 
phase difference ±i 8k z between the components in the two modes relative to this mean. 

The inhomogeneities are introduce d th rough a fluctuating part, 8K a P, of the dielectric tensor. On making the 
parabolic approximation to equation (2.1), one obtains the following equations for the propagation of the right- and 
left-hand polarized wavefields 



3 



2ik^- + Vi + ^5K++) u+ + ^SK+-U-e- 2tSkz = 0, 
oz c l ) & 

2ik— + Vi + -jSK— J u- + -^SK-+u+e 2l5kz = 0, (2.9) 



where the perturbation terms involve 



and 



8K aa , = e* a e a ,SK a P, (2.10) 



A a = e ikz e a u a (z, r), u a (z, r) = e a u a (z, r)e la5kz . (2.11) 

<T = ± 

The terms involving 5K^ and 5K |_ are zero if the inhomogcneities do not affect the polarization of the natural 

modes, that is if 5K a ^ — ctj^SKi is proportional to (K al3 ) — cr" 13 (Ki). This is the case to an excellent approximation 
if the fluctuations do not involve the direction of the magnetic field, and even for fluctuations that affect the direction 
of the magnetic field it is an excellent approximation provided that the modes are nearly circularly polarized. The 
extreme conditions under which these terms might be non-negligible are ignored here. 

Neglecting the terms involving SK^ and 5K |_, the following relations describe the propagation of the wave 

amplitude through a weakly anisotropic inhomogeneous medium 



n , 9 

2ik— - 
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u+ = 0, 

_ = 0. (2.12) 
Equations ( p. 12 ) are used as the basis for the theory developed below. 



III. SECOND ORDER MOMENTS OF THE WAVEFIELD 



In this section we calculate the average visibilities in each of the four Stokes parameters, thus determining the 
properties of the average image of a scattered source. 



A. Ensemble average visibilities 



Consider a two-element interferometer with receivers located at positions ri and r2, that measures the left- and 
right-hand circularly polarized components of the electric field. We define the visibility in a given polarization using 
the generalized second-order moment of the electric field 

7 CTCT /(z;ri,r 2 ) = (u a (z;r 1 )u* a ,(z;r 2 )) . (3.1) 

The Stokes parameters are defined in terms of the left- and right-hand circularly polarized components of the electric 

field for ri = T2'- I++ = u + u* + , I = u^u*_ 1 7_| = u + u*_ and I |_ = u-u* + . These are related to the conventional 

Stokes parameters I, Q, U and V by 

/=|(J ++ +/__), Q= +/_+), E/ = i±(J+_ -/_+), V = |(/ ++ -/__), 

I++ = I + V, J__=J-V, I+- = Q-iU, I-+ = Q + iU. (3.2) 

In this notation the Stokes visibilities are 

7(z;ri,r 2 ) = \[l++{z\ n, r 2 ) + j—{z; ri, r 2 )], 
Q(z;ri,r 2 ) = ±[ 7+ _ (z; n, r 2 ) + j-+(z; ri, r 2 )], 
U{z;ri,r 2 ) = [z] v x , r 2 ) - j- + (z; r ly r 2 )], 

V{z;r 1 ,r 2 ) = \[l++{z\ r x , r 2 ) - 7 — (z; r x , r 2 )]. (3.3) 
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Using equations ( |2.12 ) and their complex conjugates, the four visibilities obey the four propagation equations obtained 
by setting a — ±1 and er' = ±1 in 



az 



+ k 2 8K^{vi) - k 2 6K a > a ,{r 2 ) 



^aa'(z- 1 ri,r 2 )e 



iSkz(a' —(j) 



0. 



(3.4) 



Due to both the properties of the radiation from the source and the stochastic nature of the phase screen, it is 
desirable to compute the ensemble average visibility, denoted by T aa i . The ensemble average is considered an average 
over both time and over the phase fluctuations. It is normally assumed that the phase screen itself is static, with 
any perceived time- variability due to relative motion between the screen and the source-observer line of sight at some 
velocity v. This is known as the frozen screen approximation. An observer measures the visibility function at time t 
to be 7C-C-/ (z; ri + vt, r 2 + vt). However, if the system is homogeneous this visibility depends only on ri + \rt — r 2 — vt 
and is therefore independent of t. Thus the average over time is trivial and the visibility is a function of relative 
receiver separation only. 

On replacing the independent variables ri and r 2 by r = ri — r 2 , s = \{v\ + r 2 ), it follows that in a statistically 
homogene ous medium the average over the fluctuations is a function of r only. The average over the propagation 
equation (|3.4|) then leads to a propagation equation for, T a(y >(z;r) — {^ aa ' (z; r + r', r')): 



2ik 



d_ 

dz 



V r . s + fc 2 ^(r) 



iSkz(<j' — a) 



o. 



(3.5) 
(3.6) 



with V r . s = 2(d 2 /dr x ds x — d 2 / dr y ds y ), r = r x x + r y y, s = s x x + s y y. Suppose that the phase inhomogeneities are 
located on a thin screen of thickness Az. Then to first order in Az, the visibility measured a distance z from the 
screen is independent of the screen thickness. It is related to the incident visibility T aa i (0;r') according to 



k 
2nz 



T aa ,(z;r) - 

A<p aa ,{z;r') = k 2 f dz'U'iz';?') 
Jo 



d 2 r' d 2 sT r7r7/ (0;r') exp 



ik(r — r') 



+ iA(j) aa ,(r') 



(3.7) 
(3.8) 



B. The average over phase fluctuations 



The average over the random fluctuations on the screen is performed under the assumption that the fluctuations 
in both the isotropic and anisotropic terms are gaussian. The average (exp[ife]) = exp[— ((5x) 2 ) /2] applies for any 
gaussian random variable Sx. Writing 5<p a (r) = J„ z dz'(ui/c)5K crcr (z';r), we make use of this average in defining the 
generalized phase structure function: 

D aa ,{r) = 2[C CTCT ,(0) - C aa ,{r)], (3.9) 
C CTCT (r) - (S^(r')d^(r' + r)). (3.10) 

On specializing to the case where the na tur al wave modes are circularly polarized, the contribution of the isotropic 
and anisotropic fluctuations in equation (|j]) may be made explicit by introducing the notation SK aa = SKj + aKy, 
where Ki is the isotropic component of the tensor and Ky is the anisotropic component. (The cllipticity of the modes 
is determined by Kq/Kv, which is set to zero in assuming that the polarizations are circular.) The phase fluctuations 
may be separated in an identical manner: 6<p a (r ) — S(f)j(r) + a8(j)yir) with <j>i and <j>y denoting the isotropic and 



anisotropic components, respectively. Equations (3.9) are expanded as follows 

D±±(t) = Djj(r) ± 2D IV (r) + D vv (r) (3.11) 

D± T (r) = D H {v) T D IV (r) ± %(r) - D vv {v) + AC vv (0), (3.12) 

with 

Cxy = (%(r)%(r' + r)), X,Y = [I,V]. (3.13) 

The structure function Du(r) represents the effect of the isotropic phase fluctuations. Fluctuations purely in the 
rotation measure are characterized by Dyy which we call the "rotation measure structure function" . The terms Djy 
and Dyi represent the cross-correlations between in the isotropic and anisotropic phase fluctuations. 
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Perfor ming t he average over the phase fluctuations and denoting the mean anisotropic phase Skz/2 introduced in 
equation ( 2.11 ) by —<ftv the- ensemble-average mutual coherence is 



T a „i{z\T) 



(— 

\2irz 



d 2 r'd 2 s'r CTCT /(0;r') exp 



ifc(r-r')-s' D aa ,(r') 



+ i(a 



(3.14) 



Equation (3.3) is inverted at z = to obtain initial values of T aa i in terms of the Stokes parameters at the screen 
z = 0, 7(0; r), Q(0;r), [7(0; r) and y(0;r). The ensemble-average visibilities may therefore be expressed in terms of 
the initial polarization and the generalized phase structure function as follows 



(I)(z;r) 

<<9>(*;r) 

(U)(z; r) 
(V)(z;r) 



J(0;r) + V(0;r) . , 7(0; r) - 1/(0; r) . . . 

o exp[-|j> ++ (r)] + ' ; exp[-l/J._(r)], 



Q(0;r)-iE/(0;r) 2i 



2 " e 1 ' M ' 1 2 

»Q(0;r) + 77(0;r) ^ 
2 

7(0; r) + U(0;r) 



exp[-i J D+_(r)] + 



exp[-±L> + _(r)] 



Q(0;r) + iT7(0;r) 

2 ' 

iQ(0;r)-77(0;r) 



exp[-§D_+(r)], 



exp[-±£>_ + (r)], 



, J(0;r)- V(0;r) x 
exp[-±D++(r)] exp[-±D__(r)]. 



2 rL 2 /J 2 

The assumption that the statistics of the phase fluctuations are stationary implies (6^i(r)d<pv{^ + r')) 
)(j)v( r )S(f>i(r + r')), and hence = 7?y/- The visibilities reduce to 



(3.15) 



(77)(z;r) 
V(F)(z;r) 



/ e- I?v - v 'W/ 2 [(7)(0;r)cosh7? / y(r) - (F)(0;r) sinh£» /v (r)] \ 
e Dvv(r)/2-2Cvv(0){{Q/( 0;r ) COS 20y + (77) (0; r) sin20 y } 
e flvv(r)/2-2C vv (o) {_ (Q) (o ; r ) sin 20 v + (77) (0; r) cos 2^ v } 



-D vv (r)/2\ 



(3.16) 



(7)(0;r) sinh7^ /v (r) + (F)(0;r)coshJ5 /v (r)] / 



The mean values of the Sto kes p arameters I and V are equal to their incident va lues ; this may be seen by setting 
r equal to zero in equation ( |3.15 ), and noting D aa i(0) = according to equation (3J3). However, even if the initial 
circular polarization is zero, its visibility (V)(z; r) is nonetheless non-zero by virtue of the difference between 7? ++ (r) 

and D (r) at nonzero r. An interferometer may, as opposed to a single dish, therefore detect circular polarization 

from a source even if its radiation is intrinsically unpolarized. This effect was identified by Kukushkin & OPyak 



The depolarization of the linear Stokes visibilities due to stochastic Faraday rotation is manifested through the 
termexp[-2C yy (0)] §. 

In a simple model for rotation measure fluctuations in a homogeneous magnetic field (see § IVC[) , the resulting 
circularly polarized visibility (V)(z;r), is of order a times smaller than the mean intensity for initially unpolarized 
radiation. Effects of order a are far too small to be of interest for the ISM. Nevertheless, it is of formal interest to 
interpret their origin. In a medium with a homogeneous magnetic field one has 5<pv — ot<j)i, and the phase structure 
function for the right-hand polarized wavefront is (1 + a) 2 Du(r) and that for the left-hand polarized wavefront is 
(1 — a) 2 Djj(r). Thus the scale-length over which each wavefront experiences a root-mean-square phase difference of 
one radian differs slightly. This is interpreted in terms of each sense of circular polarization corresponding to different 
parameters, rdig+ and rdiff_, say, for the diffractive scales. If a is positive, one has rdig+ < raifr-j and the left-hand 
polarized visibilities extend to larger baselines than the right-hand polarized visibilities. 



IV. FOURTH ORDER MOMENTS OF THE WAVEFIELD 

The discussion in §3 refers to the ensemble averages of the Stokes visibilities. In this section we derive the variance 
of the visibilities and hence of the Stokes parameters themselves. The underlying idea is that (V 2 ) ^ and (V) = 
can lead to observable circular polarization provided that the timescale for the fluctuations in V is long compared 
with the timescale for an observation. 



A. Solution for statistically homogeneous fluctuations 



The fourth order moment of the wavefield describes the correlations of the Stokes visibilities. Using the definitions 
of the Stokes parameters in equation (3.2) the auto- and cross-correlations of the Stokes visibilities may be expressed 
in terms of the following generalized fourth-order moment of the electric field 



G 



la^'y^z; ri, r 2 , r\, r' 2 ) = (u ai (z; ri)ii CT2 (z; r 2 )<, (z; ri)i£, (z; r 2 )). 



(4.1) 



This moment describes the cross-correlation in the electric field between four receivers at positions ri, r' 1; r 2 and r 2 , 
each receiver measuring either the right- or left-hand circularly polarized component of the radiation according to the 
sign of the su bscri pt <7i, cr-2, a^, a 2 = ±1. 

Equations (2.12) are used to derive the following equation for generalized fourth-order moment (actually 16 equations 
for the 16 moments): 



2ifc^ + V 2 + V 2 - Vf - V 2 2 +ikG' aiaWi< {z-,v l ,v 2 ,T , 1 ,v' 2 ) S j 7(TlCT2CTK (z;r 1 ,r 2 ,r , 1 ,r 2 ) = 0, 



(4.2) 



with G' , a , being the ensemble average over the phase fluctuations. For o\ + cr 2 — <j[ — o' 2 7^ the following result 
holds: 

G'^, a , (z; rij r 2 , r[, r' 2 ) = D'^, (n -r[) + D'^, (z; n - r 2 ) + (z; r 2 - r[) + D'^, (z; r 2 - r 2 ) 

-^ iff2 (z; n - r 2 ) - (r' x - r' 2 ) + 2z0' y (^i + ^2 - ^ - <r 2 ), (4.3) 

where the primes on G CT1(72a ^ CT ' (z; ri, r 2 , r' 1; r 2 ), -D(z; r) and 0y denote derivatives with respect to z and the dependence 
of <fr'y on z is implicit. 



Changing co-ordinates to 



R= J(ri+r 2 +r'i+r 2 ), 
r = ri + r 2 - r[ - r 2 , 
Pi = t;0i -r 2 +r; -r 2 ), 

P 2 = ^Oi - r 2 -r'j +r 2 ), 



(4.4) 



it is evident that r ffl(T , a j^ docs not depend on R in a homogeneous medium, and this also enables us to eliminate r 
as a parameter. In view of this simplification it is convenient to change notation, replacing 7 . ia . 2Cr ' ff / (z;R, r, p l5 /t> 2 ) 
by r CTl<T2(7 / . 2 (z;p 1 ,p 2 ). Equation ( fT^ ) becomes 



2^ + 2 V Pl ■ Vp 2 + ikG'„ ia2<a , 2 (z; Pl ,p 2 ] 



r CTlO- 3 CTja- 2 ( Z i Pi; P 2 ) — 0, 



(4.5) 



G'^^ (z; Pi, P 2 ) - r>; iffi (*; P 2 ) + 2^ (*; Pi) + D'^ (z; Pl ) + D' a2< (z; p 2 ) 

-D' aia2 (z; Pl + p 2 ) - D' <a , 2 (z; Pl - p 2 ) + 2*^(cn + <r 2 - - cr 2 ). (4.6) 

The generalized fourth-order moment, r,,,^^,,' (z;p 1 ,p 2 )) describes the correlation between four receivers arranged 
in a parallelogram whose axes are Pl and p 2 (see figure 1). In an isotropic medium the co-ordinates Pl and p 2 are 
interchangeable, T aia2<y i a i (z; Pl , p 2 ) = ^a 1 a 2 a-' 1 a-' ( z : P21 Pi)i (§20-13 of p0[), but this is not the case in general in an 
anisotropic medium. 

It is useful to relate the generalized fourth-order moments to auto- and cross-correlations of the Stokes visibilities. 
The 16 terms in T ail72(7 i a i 2 are separated into four (a± — a[ = a, cr 2 = a' 2 = a' with a = ±1,<t' = ±1) that involve 
only / and V: 

T at7 > at7 >(z, Pl ,p 2 ) = (I 2 )(z, Pl ,p 2 ) + (a + a'){IV)(z, Pl ,p 2 ) + aa'{V 2 )(z, Pl ,p 2 ), (4.7) 

and four [a\ = cr 2 = a, a[ = a 2 = 0' and <x\ = cr 2 = a, a[ = a' 2 = a' with a 7^ a') that involve only Q and U: 

r***><r>(z,p 1 ,P2) = (Q 2 )(z,Pi,P 2 ) - 2ia(QU)(z,p 1 ,p 2 ) - (t/ 2 )(z,p 1 ,p 2 ), 

r a ^^ a (z, Pl ,p 2 ) = (Q 2 )(z, Pl , p 2 ) + (U 2 ){z, Pl , p 2 ). (4.8) 

The other eight terms involving cross correlations between /, V and Q, U are not discussed here. 
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B. Solution for a thin screen 



A standard approximation in scintillation theory is to assume that the phase fluctuations in the medium occur 
on a thin screen located a distance z from the observer ]30| , pT[ |. We assume that the incident wavefront is planar, 
corresponding to a point source at z = — oo so that (XF)(0, p 1 , p 2 ) is independent of the transverse spatial co- 
ordinates pi and p 2 . This implies that the fourth-order moments incident on the screen are not functions of p l or 
p 2 , and we henceforth write (XY}(0) for the incident value. 

If the source is located at a finite distance it is possible to correct for the spherical nature of the wavefront by 
making the substitutions |31[ 



ziz 2 

z\+z 2 ' 
z\ 

r. 

z\ + z 2 



(4.9) 
(4.10) 



where z\ is the distance from the source to the scattering screen and z 2 the distance from the screen to the observer. 
In particular, equation ( i.K ) implies that the length scale of fluctuations on the observer's screen is larger by a factor 
(zi + z 2 ) I z\ compared to th e pl anar case. 

The solution of equation (4.5) [j30| for a planar wavefront incident upon a thin screen located a distance z from the 
observer is 



^a 1 a 2 a' 1 a' 2 (z,r 1 ,r 2 ) 



k 

2ttz 



~ik(rx - r[) ■ (r 2 - r' 2 ) 1 



x exp 



2 L ^ r a- t a 2 cr' 1 a- 2 



(4.11) 



With this solution, we use equations ( |4.7| ) and (4.8) to assemble solutions for the correlations of the Stokes visibilities 
as follows: 



(XY)(*,n,r a ) = 




y 2nzJ J [ z 

Wn 2W IV Wvv \ ( (I 2 )(0) \ 

W IV Wn-Ww W IV (IV 2 ){0) 

W vv 2W IV W u ) V (^ 2 )(0) ) 

(W QQ -2W qu W vu \ ( (Q 2 )(0) 

W QU Wqq-Wuu -W qu (QU 2 ){0) 

\Wuu 2W QU W QQ J \ (U 2 ){0) 



(4.12) 



(4.13) 



where the arguments (r^r'j) are omitted, and we introduce 




/ W QQ \ a 

[w QU = e - 

\WuuJ 1 



e~ a cosh b + e a ~ c 

~e~ a sinh6 
e~ a cosh b — e° _c 

e' c '+d(f + f) 

id(f-f')/2 
e- c ' - d(f + f) 



a = G vv (r' x ,r' 2 ), 



b = 2G IV (r\,r' 2 ), c = 2D vv (r' 2 ), J = 2D vv (r[), 
d = D vv (t' x +r' 2 ) + D vv {t' x - r' 2 ) - 8C V v(0), 
f = explD/vW + r> 2 ) - D IV (r[ - v' 2 )]e-^ v /2, 
/' = cxp[-^ /v (r' 1 + v'. 2 ) + D IV (r[ - v' 2 )]e^ v /2, 



(4.14) 



G XY {r U T 2 ) = - {2D xy {tx) + 2D X Y(r 2 ) - D XY (ri + r 2 ) - D XY {r x - r 2 ) 



X,Y = (I,V). 



(4.15) 
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In the absence of a magnetic field, an obvious, although important point is that all the Stokes parameters scintillate 
like the total intensity. Referring back to the definitions of Djy and Dyy in §3, the absence of Faraday rotation terms 
implies that Djy, Dyy and cj)y are zero, leaving only Wu, Wqq and Wuu nonzero. In particular, one has Wu = 1 
and Wqq + Wuu — 1> so the correlation functions differ only by a multiplicative constant ((XY)(Q)): 



{XY)(z, ri ,T 2 ) 
XY 



d 2 r' 1 d 2 r2 exp 
[II,IV,VV,QQ,UU,QU]. 



ik(ri -r[) ■ (r 2 - r 2 ) 



G//(ri,r 2 ) 



(4.16) 



C. A simple model for rotation measure fluctuations 



A simple model is when the fluctuations occur only in the electron density, with the magnetic field being uniform. In 
this case the structure functions Dyy and Djy are related to Du by the parameter a — Q e /ui, where f2 e is the electron 
cyclotron frequency and uj is the angular frequency of the radiation. In particular, one has Dyy /a 2 = D\y ja = Du. 
Taking a typical value of the magnetic field in the ISM of 3 /xG and an observing frequency oi v = 1 GHz, a is of order 
10 \ 

In this case we expand equations (4.14) in terms of a, where Djy /Du is of order a and Dyy / Du is of order a 2 . 



Retaining terms to second order in a the second order correlations are given by equations ( 4. 12j ) with 




l + 6 2 /4-c/2 
6/2 

6 2 /4-a + c/2 
1 + a-c' + „ e -8CVv(o) 



9£ 

/j e -8CW(0) 



1 + a 



5e 



-8CVv(0) 



(4.17) 



g = cos 4(f>y + i sin 4(j)y (Z3/y(r' 1 - r' 2 ) — D IV (r[ + r' 2 )) 

+ l - cos40 y { [D IV (r[ - r' 2 ) - D IV (r[ +v' 2 )f 

-Dyy^ + r' 2 ) - Dyy^ - v' 2 ) - GW(ri,4)} 

h = sin 40y — icos40y {Djy^ — r' 2 ) — Di V (r\ +r' 2 )) 

sin4(/)y r 2 
— 1 [D IV (r 1 - r 2 ) - D IV (r 1 + r 2 )] 

+D vv (r , 1 + v' 2 ) + D vv {v\ - v' 2 ) + GW(riy 2 )}. 



+ - 



(4.18) 



Equations ( 4.17 ) and ( 4.18 ) are used in the discussion below. 



V. IDENTIFICATION OF EFFECTS 



In this section we discuss the interpretation of the equations ( 4.1 2j ) and (4.17), concentrating on the fluctuations in 
circular polarization and the total intensity. 



A. Circular Polarization 



We are concerned with the creation of a circularly polarized component and so we assume the source has no intrinsic 
circular polarization, which corresponds to assuming (V 2 )(0) = and (/V)(0) = in equation (4.12). The variances 
in the visibility and in the Stokes V visibility are then 
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(I 2 ){z,r x ,r 2 ) = <7 2 }(0) (J^j J dMr' 2 exp 

x{l + \G IV (^' 2 f 









( cPy^t^ exp 





' ikjrx - r[) ■ (r 2 - r 2 ) 

^ II V r l> r 2) 

-D vv (r' 2 )}, 

ikjrx - r[) • (r 2 - r 2 ) , , 
uii(r 1 ,r 2 ) 



(5.1) 
(5.2) 



x {-Gw(ri,r 2 ) + D vv (r' 2 ) + |Gjv(ri, r 2 )| 2 } . (5.2) 

For convenience in the discussion below, we label the contributions to (V 2 ) due to Gyv, Dyv and |G/y| 2 as (V 2 )\, 
(V 2 )2 and (V 2 ) 3 in (5.2) respectively. 



1. Variations in the total intensity 



Before discussing the interpretation of equation ( |5.2| ), we review the subject of intensity variations in an isotropic 
medium because many of the approximations and definitions are relevant to the anisotropic case. In the strong 
scattering regime, intensity variations are due to diffractive scintillation, caused by interference between subimages 
over the scattering disk, and refractive scintillation, due to refractive focussing and defocussing of the entire scattering 

disk i-. 

Intensity variations in an isotropic medium correspond to the first term in the curly brackets in equation ( p.l| ). 
Following 133] this term may be written in the form 



(P)u(z;t u t 2 ) = 



= /j L V 



\2-kz J 



d 2 r' 1 d 2 r 2 (7 2 )(0) exp 



ik(ri 



(r 2 



[e~ D ^ [1 + Q(r 2 , r[) + ■■■] + e~ D W [1 + n(r[, r> 2 ) + •••]}, 



n(r 1( r 2 ) = Dnin + r 2 )/2 + £> /7 (n - r 2 )/2 - D„(ri), 



(5.3) 
(5.4) 



where the subscri pt 7 7 signifies that only isotropic phase fluctuations are retained. The approximation made in 
deriving equation (5.3) follows by recognizing that, for strong scattering, the intensity fluctuations are dominated by 
r egio ns where either r[ or r 2 are small and then assuming r 2 3> r% or n ^> r 2 for these regions respectively. Equation 
(5.3) is written compactly in the form 



(7 2 )„(z,n,r 2 ) = [2 + 2e(z,r 1 ,r 2 )](7 2 )(0). 



(5.5) 



The term £ is due to the contribution of f2, and is associated with large scale focussing and defocussing of the entire 
scattering disk, giving rise to so-called "refractive scintillation" H. The explicit expression for £ is 



f d 2 q 

£0;ri,r 2 )=2 J ^-p-$(q){l-cos(q-p 1 - g 2 r|)}exp[iq-p 2 -7> // (p 1 -qr|)]. 



(5.6) 



This term is less than unity in the strong scattering regime ||. The power spectrum of phase inhomogeneities, < I ) (q), 
is defined by 



D H (r) = 

A power law spectrum ^(q) has the form 

*(q) - Qii<T l \ 



rf 2 q 



$(q)[l 



9min < q <q n 



(5.7) 



(5.8) 



where q m i n and q m & x correspond to the inverses of the outer and inner scales of the scattering medium, r out and ri n , 
respectively, and Qn is a constant. The familiar case of Kolmogorov turbulence corresponds to (3 — 11/3. Relating 
the definition of the phase structure function Djiijc) provided by |33| to the power spectrum of phase inhomogeneities 
p3, one has 



Qu = 8n 3 ALC%r 2 e \ 



(5.9) 
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is the classical radius of the electron, AL is the path length through the scattering medium and Cfr is the 



where r, 

electron density structure con stant. The form ( |5.g| ) for ^(q) is used extensively below. 
Inspection of equation (5.5) reveals that the variance of the intensity scintillations, 



(I 2 ) (z, 0, 0) - (I 2 ) (0) = [1 + 2£(z, 0, 0)] (I 2 ) (0) 



(5.10) 



contains contributions from two terms. The term 2£(z, 0, 0) (I 2 ) (0) is associated with refractive scintillation, while the 
term of magnitude (I 2 ) (0) is identified with diffractive scintillation. 



2. Circular polarization correction terms 



We now consider the three terms in equation (^j^) , assuming that fluctuations in the rotation measure are determined 
by fluctuations in the density alone, implying Du = Djy/a — Dyy jo? with a = Q, e /uo. 

The contribution from the first term, denoted (V 2 )i, comes from the re gion s wher e ei ther or r' 2 are small. Since 
Gvv — o> 2 Gu, this term may be approximated by referring to equations (5.3) and (5.6): 



i( k 
— ce - — 

\2-kz 



dr' x dr' 2 {I 2 )(ti) exp 



ik{r x - ri) ■ (r 2 - r' 2 ] 



{exp[-D // (ri)]n(r^r , 1 )+exp[-£) // (r^)]n(r' 1) r^)}, 



which may be rewritten as 



(y 2 ) 1 (^r 1) r 2 )=a 2 (/ 2 )(0)[^(z ) r 1 ,r 2 )+^,r 2) r 1 )], 



(5.11) 



(5.12) 



with £ given by equation (|5.6|). Production of circular polarization due to the term (V 2 )\ is associated with the 
correction term in equation Q5.6| ) due to refractive intensity variations 
It is possible to de rive an expression for (V 2 )i from equation (5.11 



brackets in equation (5.11). The exponential term is only large for sma 
for small ri/r 2 : 



) directly. Consider the first term in curly 
1 r[, so it is appropriate to expand f2(r 2 , r^) 



(— 



0- 2 /_/\2 

(/3-3)(/3-2) U 



(5.13) 



Specializing to the case ri = r 2 = 0, the terms involving r[ and r 2 in curly brackets in equation (A3) are interchange- 
able. One then has 



(^ 2 ) 1 (z,0,0) = 2a 2 (/ 2 )(0) 



k 
2?rz 



dr' 1 (/?-3)(/3-2)r 2 



(-5- 



13-2 



dr' 2 exp 



ikr[ ■ r' 2 



//3-4 



We evaluate the integral over r 2 using K 



J d 2 r' 2 exp [ix • r 2 ] r\ = -KV+ 2 



£(7/2 + 1) 

r(-7/2) • 

which is valid for 7 > — 2. The remaining integral over r[ yields 



-7-2 



(F 2 )!(z,0, 0) = a 2 (/ 2 )(0)2^- 3 7r- 1 (/3 - 3)09 - 2) 



^diff 



2/3- 



r(2-/3/2) V/3-2 



(5.14) 
(5.15) 
(5.16) 



Since the mod ulation index due to refractive intensity variations, m re f, is of order (rp/rdiff) 4_/3 (see, e.g 
equation ( ^3 .16 ) implies a degree of circular polarization of order am rc { (J 2 ) 1 / 2 . Appendix A details the contribution 
of (V 2 )i when the rotation measure fluctuations are not necessarily proportional to the isotropic phase fluctuations. 

The contribution due to the term (V 2 ) 2 is also evaluated by representing the rotation measure structure function 
in terms of its power spectrum. Evaluating the integrals over and r' 2 in equation (5.2), the variance in circular 
polarization due to the term Dyvi^) is 
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(V 2 ) 2 (z,r 1 ,r 2 ) = 2a 2 (I 2 )(0) / d 2 q 



1 



g(q) 



exp[-£>(n)] 



- {e 4q r2 exp[- J D(r 1 - r|q) + e - 4q r2 exp[-£»( ri + r|q)} 



(5.17) 



For arbitrary ri and r 2 this expression is, in general, complex, and unlike (V 2 )i(z, r ±, r 2 ) , it is not symmetric under 
interchange of ri and r 2 . The circular polarization is a real quantity, and equation (5.17) can be directly expressed 
in terms of the Stokes parameter V only for r 2 = 0, in which case the expression is real, as required. 

Th e contribution from (V 2 } 2 to the variance in circular polarization is obtained by setting ri = r 2 = in equation 
( 5.17 ) . The integral is approximated using 



[1 - exphC(qri)]] « ( f fr K i < < lArCf ' 

L L V t>U1 {1, 1/rref < q < ?max, 



to yield 



/8-2 



In 



ymin'ref 



(5.18) 



(5.19) 



where we assume that the outer scale is much larger than the refractive scale (Lq 3> r re f). Rewriting Qn in terms of 
rdiff, one finds the magnitude of the circular polarization: 



(V 2 ) 2 (z,0,0) = a 2 (I 2 )(0)f((3) 



( rp 
V^diff 



2/3-4 



1 



0-2 



In 



1 



'Zmin? rcf 



(5.20) 



with 



= pg-i r(/3/2) 
r(i-/3/2)' 



(5.21) 



for /3 < 4. For /3 = 11/3 one has f{0) » 0.89. The contribution of the circular polarization due to this term scales as 
a 2 (r^/r dss )P- 2 . 

The contribution of (V 2 }^ in a homogeneous magnetic field is evaluated by appealing to the arguments presented 
in the evaluation of {V 2 )\. For strong scattering, only the regions near r[ « and r' 2 w contribute. Expanding for 
r[ <C r 2 and r 2 <C one has 



(V%{z 1 T 1 ,v 2 ) = a 2 (I 2 ){0)[ — 



d 2 v' 1 d 2 r' 2 exp 



ifc(rx - ri) • (r 2 - r 2 ) 



x{e- D ^n a (^ as ri) + e-^)n a (riy 2 )} 



(5.22) 



Consider the second term in curly brackets in (5.22). The exponential ensures that this term is only large for small 
r' 2 , so one may expand n 2 (r' 1 ,r 2 ), defined by equation (|5.4[), for small r 2 . To second order in r' 2 /r[ one has 



2(/3-2) 



(/3-2) 2 (/3-3) 2 r 2 4 
4r< 4 



(5.23) 



Recognizing that the two terms inside the curly brackets in equation (5.22) are identical under the replacement 
r i *~* Y 'ii provided ri = r 2 = 0, one then has 



(V' 2 ) 3 (»,0,0) 



(— 

\2-kz 



d 2 v\d 2 Y 2 exp 



ikr[ ■ r' 2 



-D n (ri) 



rfl? 2 (r 2 )(/?-2) 2 (/3-3) 2 



Now we may evaluate the integral over r' 2 using equation ( 5.15 ) . Then equation ( 5.22 ) reduces to 



(5.24) 



(5.25) 



Now, using 
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J^dnrJeM-ir/r^f- 2 } =r^-^T (j^j , (5.26) 
which is valid for (3 > 2, one has 

<U 2 ) 3 = a 2 (I 2 )(0)((3 - 2)09 - 3) 2 2 2 ^ 7 [^Y^ T (^f) f^Jy (5-27) 

For strong scattering, the term (rp /^diff ) 4 ' 3-16 is less than unity. In particular, for a Kolmogorov spectrum of turbulent 
fluctuations (/? = 11/3) the root mean square degree of circular polarization due to this term scales as a^dig/rp) 2 / 3 , 
less than unity for strong scattering. 

The dominant term is (V 2 )2 S> (V 2 )i, (V 2 )3 as this is the only term that allows -J (V 2 ) / (I) a to be much greater 
than unity. We retain only this dominant term in the discussion below. 



3. Cross Correlation 



So far we discuss only the variance (V 2 ) in Stokes V. The correlation function (IV) contains additional information 
about the propagation-induced circular polarization. 

The form of the expression for (IV)(z, rx, r 2 ) in equation fl4.12 ), in th e abs ence of any intrinsic circular polarization, 
is of the same form as the expression for (V 2 )i, yielding (cf. equation (5.11)) 



(IV)(z,0,0) = (l 2 )(z,0,0) 



k 
2-kz 



d 2 r' 1 d 2 r' 2 exp 



ikr[ ■ r' 2 



x {exp[-£(r;)]f! /y (r' 2 ,ri) + exp[-D(r' 2 )]Q IV (r' 1 ,r' 2 )} , 
J y(r 1 ,r 2 ) = D IV (r 1 + r 2 ) + D IV (r 1 - r 2 ) - D IV (ri). 

For a homogeneous magnetic field, one has Q.iv(yi,Y2) = a;fi(ri,r2) and the cross correlation is 

</U)(z,0,0) = 2a£(z,0,0)(/ 2 )(0). 



(5.28) 
(5.29) 



(5.30) 



Since 0, 0) is of order unity, it follows that the cross correlation is of order a, rendering this correlation too small 
to be observed. 



VI. DISCUSSION AND CONCLUSIONS 



Scintillations in a magnetized ISM necessarily lead to a small degree of circular polarization (CP) even for an 
unpolarized source. This occurs because the refractive indices for the two natural wave modes are slightly different, 
leading to a variety of small effects that are different for the two opposite CP components. Under the assumption 
that the fluctuations involve only the plasma density these effects are characterized by a single small parameter, 
a = Q e /u;, which is the ratio of the electron cyclotron frequency at the scattering screen to the wave frequency. For 
typical parameters in the ISM, a ~ 10~ 8 is too small for the induced CP to be of practical interest unless the effect 
can be enhanced in some way. 

The effect that we identify as of possible practical interest appears in the variance of the Stokes parameter V, 
denoted (V 2 ). Provided that the time scale of observation is short compared with the time scale on which (V 2 ) 
changes, a net circular polarization of order (V 2 ) 1 / 2 / (I) should be observed. In evaluating (V 2 ) we separate it into 
three terms, (V 2 )i,i = 1,2,3, and find that each of them is of the form (V 2 )i = a 2 (I 2 )(0)Ai(rp / Vdiff) ai , with the Ai 
all of order unity. For strong scattering, rp ^> rviiff implies (V 2 ) 3> a 2 (I 2 ) provided that <Zj is not too small. For 
the term (V 2 )2, with a-i = 2/3 — 4 = 10/3 for a Kolmogorov spectrum (3 = 11/3, the numerical factor is large when 
the scattering is strong. This term arises from scattering due to the largest structures in the postulated power-law 
spectrum of fluctuations. 

An obvious prediction of this theory is that the observed degree of CP should reverse randomly on the time 
scale over which (V 2 ) changes. However, the simplifying assumptions we make need to be reconsidered in making a 
realistic estimate of both the magnitude and time scale of the fluctuations in V. In particular, we assume that both 
the fluctuations that lead to scintillations in the usual sense and the fluctuations in rotation measure that cause a 
separation in the right and left polarizations occur at a single phase screen. The assumption that these two effects 
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occur together is unnecessarily restrictive. The effect that we describe requires (a) that the wavefront be rippled, and 
(b) that the right and left polarizations be separated, but there is no need for these two effects to be due to the same 
structures in the ISM. Indeed, the dominant contribution to (V 2 ) is due to the large-scale structures in the assumed 
turbulence, and in fact any large-scale magnetized structure can act as a "Faraday wedge" in refracting the right and 
left polarizations into slightly different directions. In principle, the separation between the two circularly polarized 
rays can become arbitrarily large at arbitrarily large distances from such a Faraday wedge, leading to non-overlapping 
right and left hand polarized images, and so to (V 2 ) ~ ((7 — (I)) 2 )- We propose to develop this idea further and 
apply it to the interpretation of observed CP elsewhere. 

Furthermore, an exact theory requires knowledge of the turbulent fluctuations in both the orientation and magnitude 
of the magnetic field (here assumed to be homogeneous), and its correlation with density fluctuations. The turbulence 
in the ISM is expected to be highly anisotropic with respect to the magnetic field lines. As remarked in 36 1, the 
fact that the VLBI images of scattcr-broadcned images are typically only slightly elliptical (with axial ratios at most 
~ 2:1) implies that the radiation propagates through many regions with different magnetic field orientations. This is 
expected to enhance the production of scintillation-induced CP. 

We conclude (1) that scintillation-induced CP must occur due to propagation through the magnetized ISM, (2) that 
the simplest estimate that it is of order a ~ 10 -8 and so would be unobservably small, is correct only for the mean 
value of (V)/(I), (3) that the variance (V 2 ) / '((I — (I)) 2 ) can be much larger than a 2 , and (4) that a more general 
model that treats the scintillations and the separation of the opposite circularly polarized rays ( "Faraday wedge" ) in 
different ways is needed in formulating a theory that will be a useful basis for comparison with observations. 
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APPENDIX A: THE CONTRIBUTION OF (V 2 )i 

We evaluate the contribution of the (V 2 )i term when the structure functions Djj and Dyv are n °t necessarily 
proportional to each other. Writing 



Dvv(r) = 

the contribution of {V 2 )\ is 

(F 2 ) 1 ( 2 ,r 1 ,r 2 )= / dr[dr' 2 (I 2 )(0)e X p 

where we write 



fdiffR 



Pv-2 



ik(r 1 - ri) • (r 2 - r 2 ) 



x {expl-DiiWnvvfayi) + expI-D^r^fW^ra)} . 
£Vv(r 2 , ri) = D vv {n + r 2 )/2 + D vv (n - r 2 )/2 - D vv (r 2 ). 



(Al) 

(A2) 
(A3) 

(A4) 



Consider the first term in curly brackets in equation (A3). The exponential term is only large for small ri, so it is 
appropriate to expand f2yy(r 2 ,ri) for small ri/r 2 : 



tt vv (r 2 ,ri) 



T2 

r&sm. 



0V-2 



(fiv-3)(fly-2) ( H 
Vr 2 



(A5) 



Specializing to the case ri — r 2 = 0, the terms involving r' ± and r' 2 in curly brackets in equation (A3) are interchange- 
able. One then has 



(V- J ),( 3 ,0,0) = 2(7 2 )(0) (JL\ J dr[(f} v - 3)tflv " 2)r? 



p v -2 



dr' 2 exp 



ikr[ ■ r' 2 



I0v -4 



(A6) 
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We evaluate the integral over r' 2 using |U| ). The remaining integral over is evaluated to yield 



(y 2 )i(2,0,0) = {/ 2 >(0)2^- 3 7r- 1 (^ - 3)(/V - 2)r|* 



rdiffR diff r(2-/? y /2) 



0-2 



(A7) 



Setting Uyy(r) = a 2 Dn(r), one then has /3y = (3 and r dl " 



2-/3 _ „,2 2-/3 



ffR 



err 



difn ' anc ^ equation (A.7) is then proportional to 



the variance due to refractive scintillation, of order (rp/rdiff) 8 2/3 which is consistent with equation (5.16). 
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FIG. 1. Positions of the receivers in calculating fourth-order moments (cf. Ishimaru 1978). 
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